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1.  GROUND  CLUTTER  PROPERTIES 

This  section  examines  some  important  properties  of  ground  clutter  that 
influence  the  design  of  a  canceler.  Measurements  and  data  are  for  the  Norman 
Doppler  radar  which  has  a  medium  to  small  clutter  problem  as  compared  to  airports 
in  the  vicinity  of  major  cities. 


1.1  Ground  Clutter  Reflectivity 

Surrounding  ground  influences  radar  observations  in  two  ways.  Returns  from 
the  ground  through  the  mainlobe  or  si  delobes  create  strong  spectral  power  around 
zero  velocity  which  contaminates  the  spectral  moment  estimates  of  weather  echoes.  •; 

Also,  tall  objects  interfering  with  the  ri  .r  beam  act  as  blockages  which  may  3 

totally  prevent  the  illumination  of  weather  beyond.  v 

<: 

Figure  1.1  illustrates  the  blockage  surrounding  the  Norman  Doppler.  Note  that  | 
there  are  not  many  significant  blockages  above  0.5°  in  elevation,  and  that  most  of  \ 

the  blockages  extend  to  less  than  0.2°.  Thus,  we  have  chosen  0°,  0.4°  and  0.8°  \ 

in  elevation  to  be  the  angles  for  data  collection.  1 

a 

Three  different  terrains  characterize  the  area  around  the  Norman  site:  \ 

urban  Oklahoma  City  to  the  no^th,  forests  to  the  east,  and  prairies  to  the  southwest.  | 

Clutter  characteristics  from  these  three  areas  are  separately  analyzed.  | 

The  map  of  the  Norman  ground  clutter  for  the  three  elevation  angles  is  shown  | 

on  Figures  1.2,  1.3,  and  1.4.  These  displays  are  in  terms  of  equivalent  reflec-  | 

tivity  factor  (dBZ)  and  only  reflectivities  above  a  threshold  (indicated  with  the  ^ 

cursor)  are  depicted.  Total  areas  of  clutter  above  a  fixed  reflectivity  for  this  ;: 

data  sec  are  plotted  on  Fig.  1.5.  f 

I 

One  can  use  Figures  1.2  through  1.5  to  roughly  assess  the  performa  ice  of  a  | 

ground  clutter  canceler,  which  is  explained  in  following  paragraphs.  Lot  us  model  | 

the  ground  clutter  and  the  weather  signal  with  Gaussian  spectra  and  corresponding  3 

autocorrelation  functions.  % 


Figure  1.1  Survey  from  Norman  Doppler. 
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where  C  is  the  clutter  power  at  a  certain  range,  a  the  width  of  jthe  clutter 
Doppler  spectrum,  and  S  is  signal  power  and  a  spectrum  width  at  the  same  range. 

Va  is  the  unambiguous  velocity,  n  the  lag  number  while  F  indicates  the  Fourier 
transform.  Even  though  the  Gaussian  model  is  a  simplification,  .it  allows  a  wide 
range  of  clutter  spectra  to  be.  easily  modeled.  Moreover,  superposition  of  Gaussian 
spectra  can  produce  more  complicated  clutter  for  which  still  the  strongest  one 
is  dominant  and  must  be  filtered. 


Throughout  this  study  we  will  assume  that  either  a  Fourier  transform  or  a  pulse 
pair  algorithm  is  used  to  estimate  the  mean  Doppler  and  the  Doppler  spectrum  width. 
If  the  Fourier  transform  is  used,  then  the  filtering  would  be  done  by  properly 
weighting  the  time  series  data  (inphase  and  quadrature  components)  prior  to  the 
transform.  For  the  pulse  pair  algorithm,  the  most  natural  filter  would  be  a 
recursive  one,  although  a  Fourier  transform  could  be  used  twice  (once  forward  to 
obtain  and  reject  the  coefficients  near  zero  and  then  back  to  regenerate  the  time 
series). 


In  Appendix  A,  it  is  shown  that  when  pulse  pair  processing  is  used  the  ratio 

Sps 

8  =  rTT 
Cpc 


where 


(1.3) 


1 


-7r2ac2/2va2 

Pc  =  e 

must  be  larger  than  10  dB  in  order  to  maintain  the  mean  velocity  error  below  , 

1  m*s-Tfor  any  velocity. 

From  Figures  1.2,  1.3,  and  1.4  we  see  that  ground  clutter  larger  than  40-45  dBZ 
is  sparce  enough  so  that  weather  signals  of  similar  magnitudes  should  be  discernible. 

If  weather  echoes  are  50-55  dBZ,  then  meaningful  estimates  of  mean  velocity  (see 
Appendix  A)  should  be  possible  and  due  to  the  spacial  continuity  of  weather  (gust 
fronts,  etc.)  observation  throughout  the  whole  PPI  could  be  made.  But  gust  fronts 
have  low  reflectivities  and  in  order  to  measure  their  velocities,  one  needs  to 
cancel  the  ground  clutter.  For  accurate  mean  velocity  estimation  via  the  pulse 
pair  processor,  one  needs  a  10  dB  ratio  of  signal  to  clutter  residue  and  15  dB 
ratio  if  spectrum  width  is  estimated  (Appendix  A).  Thus  the  canceler  must  reduce 
the  residue  at  least  10  dB  below  signal  power.  Therefore,  a  50  dB  canceler  would 
allow  moment  estimation  up  to  about  40  dB  of  C/S  ratio. 
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Figure  1.2.  Ground  clutter  map  at  0°  elevation.  Equivalent  reflectivity  factor  (dBZ) 
above  some  threshold  is  displayed.  Range  mark  is  at  40  km.  The  thresholds  are  in¬ 
dicated  with  the  cursor's  position.  Categories  of  displayed  reflectivity  in  dBZ  are 
from  9  to  over  67.  Rote  a  false  indication  of  uniform  reflectivity  (about  8  km 
in  diameter)  centered  at  the  radar,  that  decreases  in  circular  steps.  This  is  due 
to  receiver  saturation  which  creates  a  constant  power  at  the  digital ,  integrator 
input.  When  that  constant  power  is  multiplied  with  the  range  squared,  r%,  to 
obtain  reflectivity,  a  monotonic  decrease  occurs  as  r  gets  smaller.  Tnerefore,  at 
this  elevation  angle  the  radar  is  totally  blind  to  4  km  in  range.  Peak  transmitter 
power  was  200  kW. 
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Figure  1.3 .  Same  as  1.2  except  the  eleva¬ 
tion  angle  is  0.4° 3  and  the  threshold 
in  the  top  left  figure  contrary  to 
cursor  indication  is  at  -9  dBZ. 


Figure  1.4.  Same  as  1.2  but  the  elevation  is  .8°  and  the  radar  is  blind  to  about 
2  Ian. 


5 


This  margin  is  also  needed  between  the  spectrum  peak  and  its  sidelobes  if  the 
Fourier  method  of  eliminating  clutter  is  employed.  This  is  so  because  clutter 
spectrum  widths  are  typically  about  0.2  m«s  so  that  the  spectrum  peak  ratio  from 
(1.1)  and  (1 .2)  is 

Cav  __  TOC 

Soc  *  S  (1.4) 

for  a  a  of  2  m-s~\  and  we  note  that  50  dB  of  clutter  spectrum  cancellation  would 
allow  signals  with  powers  40  dB  below  clutter  powers  C  to  be  seen. 

Cumulative  areas  of  ground  clutter  obtained  for  the  three  elevation  angles 
(Figures  1.2,  1.3,  and  1.4)  are  shown  on  Figure  1.5.  We  notepthat  the  area  of 
ground  clutterpwith  reflectivity  larger  than  30  dBZ  is  400  knr  at  0°,  300  km  at 
0.4°  and  50  knr  at  0.8°.  For  the  Norman  radar  gust  front  echoes  with  -10  dBZ  equiva¬ 
lent  reflectivity  factor  are  above  noise  up  to  30  km  which  is  also  the  range 
of  ground  clutter  extent.  Therefore,  a  50  dB  canceler  would  enable  us  to  see  such 
weak  echoes  through  most  of  the  PPI  except  in  the  above  mentioned  areas. 

Measurement  of  ground  clutter  cross  section  per  unit  area  yielded  values  about 
40  dB  below  1  nr/nr  with  a  typical  example  in  Table  1.1.  Note  that  these  data  are 
about  10  dB  lower  than  the  ones  found  in  Nathanson  (1969),  which  is  not  unusual 
considering  the  variability  of  clutter  characteristics. 

Table  1.1 


Range  (km)  Average  clutter  cross  section 

(dB  below  lm  /hi  ) 

8.3  -  11.9  -39.96 

13.1  -  16.7  -51.84 

17.9  -  21.5  -40.01 

22.6  -  26.3  -68.26 


1.2  Doppler  Spectrum  Width  of  Ground  Clutter 

Two  mechanisms  produce  broadening  of  the  ground  clutter  spectrum.  These 
are  movement  of  trees  under  the  influence  of  wind  and  the  rotation  of  the  antenna 
while  data  are  being  collected.  Groginsky  and  Glover  (1980)  have  fitted  the  data 
of  Nathanson  for  the  width  due  to  wind  with  the  following  regression  line: 

%'8-  l°'3vw'-2  (m-s-’l 

where  the  wind  velocity  v, ,  is  in  m.s"1. 

w 

Antenna  rotation  at  a  rate  a  contributes  (Nathanson  1969) 

2  /a  A  cose  \ 

aca  y  2tt6-j  J  £n2 


(1.5) 


(1.6) 
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where  0,  is  the  one-way  half-power  beamwidth,  9  the  elevation  angle  and  X  the 
wavelength.  The  curves  for  o  and  cr  are  plotted  on  Figure  1.6  on  which  one 
can  observe  that  a  1°  beamwidlrl  antenna  rotating  at  3  rpm  produces  a  width  of 
.23  m*s'J.  At  wind  speeds  of  10  m*s"  (20  knots),  the  width  is  broadened  by 
.12  m*s"  .  Thus  the  total  rms  clutter  width  would  be 

'  O  O 

a  =/a  +a  oc  -1 

c  cw  *  cot  —  .25  m*  s 


Spectral  analysis  of  time  series  data  collected  with  the  Norman  radar  confirms 
this  finding.  Some  typical  spectra  (Figure-.!.?)  show  what  is  to  be  expected  at 
Norman  with  rotation  rate  1°  s"  and  10°  s  .  The  differences  in  widths  at  these 
two  rates  are  difficult  to  measure  because  noise  and  sixty  cycle  harmonics  mask 
the  weak  spectrum  skirts.  Still, we  note  that  the  third  spectral  points  both  side 
of  zero  are  higher  at  the  10°  s"  rotation  rate  due  to  more  broadening.  If  5 
spectral  coefficients  about  zero  were  eliminated  (i.e.,  -1.06  to  1.06  m-s"),  about 
50  dB  of  canceling  would  be  achieved. 


WIND  SPEED  v*  (m  s'1) 

5  10  15  20  25 


ANTENNA  ROTATION  RATE 


Figure  1.6.  Clutter  spectrum  widths  aflq  due  to  antenna  rotation 
(scale  on  lower  axis)  and  oaui  due  to  wind  speed  (upper  axis) . 
Dashed  line  is  through  two  estimated  data  points  (see  text). 
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Figure  1.7.  Plot  af  clutter  power  spectra  (128  points)  in  dB  helow  peak  versus 
velocity  in  m-s ~  using  from  top  to  bottom  rectangular,  von  Harm  and  Hamming 
windows.  The  antenna  is  rotating  at  the  rate  of  1.0 °  s'1  and  10.  0°  s  .  Spectral 
powers  are  marked  with  X  's  and  a. five  point  running  average  is  drawn  for  visual 
clarity.  The  peaks  at  ±  10  m‘8~  are  third  harmonics  of  60  Hz  picked  up  by  the 
video  channel  prior  to  the  A/D  converter. 


Figure  1.8  Bias  of  the  width  estimate 
due  to  the  finite  time  window  associ¬ 
ated  with  Fourier  spectral  processing. 
To  obtain  the  bias  in  m-s~L ,  multiply 
the  ordinate  values  with  the  Nyquist 
interval  (68  m-s  for  the  Noman 
Doppler  radar). 


Estimation  of  narrow  spectrum  widths  is  complicated  because  the  bias  effect  due 
to  finite  data  lengths  may  be  much  larger  than  the  width  to  be  measured.  For 
instance,  the  von  Hann  window  with  128  samples  contributes  on  the  average  0.3  m*s“1 
of  bias  to  the  measured  width,  while  with  1024  points  the  bias  is  .04  m-s  . 

(See  Figure  1.8.)  For  practical  reasons  histograms  of  spectrum  width  were  computed 
utilizing  a  128  point  discrete  Fourier  transform.  To  eliminate  the  effects  of 
noise  and  60  Hz  harmonics  only  ten  spectral  lines  around  zero  were  considered. 

Our  results  (Figures  1.9-1.11)  have  very  little  terrain  dependence  probably  because 
the  winds  were  light  (4.63  m*s"  ).  At  a, rotation  rate  of  10°s”  the  mean  widths 
are  about  four  times  larger  than  at  l°s  ,  rather  than  10  times  as  suggested  by 
(1.6).  We  thus  must.conclude  that  biases  are  significant.  Roughly  the  unbiased 
mean  width  at  10°-s"  would  be  about 

/A2  -  .3Z  *  .26  m-s'1  , 

which  is  plotted  on  Figure  1.6.  The  formula  (1.6)  predicts  .17  m-s”  for 
broadening  due  to  antenna  rotation  and  the  higher  measured  value  can  be  attributed 
to  wind  effects  and  other  broadening  mechanisms  not  included.  An  independent  check 
of  the  bias  in  the  measured  widths  was  made  from  the  scattergram  of  the  widths 
obtained  with  the  128  point  transform  versus  1024  point  transform. (Figure  1.12). 

To  a  .4  m-s"  width  (from  the  128  point  FFT)  corresponds  .28  m-s"  from  the  1024 
point  FFT,  which  is  quite  close  to  the  .26  m-s"  obtained  from  the  theoretical 
bias  curves. 

A  sensible  unbiasing  of  the  width  at  l°-s_1  rotation  is  difficult  because 
the  predicted  window  bias  is  often  larger  than .the  total  measured  width. 

Thus,  we  accept  the  measured  value  of  .]  m-s'  .  In  sumnary,  the.Doppler 
spectrum  width  is  between  .1  and  ,b  m-s"  with  the  mean  at  10°  s"  of  about 
.25  ms"  .  The  relationship  of  two  measured  points  to  the  theoretical  curves  is 
shown  on  Figure  1.6,  and  we  note  that  the  slope  of  the  measured  curve  is  close 
to  the  theoretical  prediction. 

Now  let  us  assume  that  L  dB  of  clutter  cancellation  (in  the  spectrum  domain) 
is  desirable  with  the  filter  having  a  passband  velocity  v  .  Then  for  an 
ideaiized  characteristic  as  on  Figure  (1.13),  we  can  findpan  approximate 

relationship  between  v_  and  the  clutter  spectrum  width  o  by  requiring  that 
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Figure  1.10.  Same  as  1.9  but  returns  are  from  the  prairies  southwest  of  the  radar 
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Figure  1.11.  Same  as  1.9  but  the  returns  are  from  the  wooded  area  around 
Lake  Thunderbird  east  of  the  radar. 
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CLUTTER  WIDTH  IN  m  »*' {m  «  1024) 


Figure  1.12.  Scattergram  of  clutter  spectrum  widths 
obtained  from  a  128  point  Fourier  transform  versus 
width  values  from  a  1024  point  transform.  Both 
sequences  were  weighted  with  the  von  Harm  data 
window.  The  solid  line  represents  a  least  squares 
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Figure  1.13.  Idealized  frequency  (velocity)  j 

response  of  the  high  pass  filter.  The  -j 

stopband  width  is  2  vs  and  passband  ;j 

cut-off  width  is  2  v  .  -i 

P  *1 

which  is  about  1  m-s-1  at  o  =  .2_m*s-1.  However,  for  widths  from  .1  to  .5  m-s”\  j 

the  v  should  span  .48  to  2.4  m>s"  .  This  indicates  that  the  clutter  filter  should  | 

have  a  notch  of  variable  width  dictated  by  the  antenna  rotation.  | 


Both  the  recursive  filter  and  the  FIR  filter  require  certain  time  to  achieve  -J 

clutter  cancellation.  For  a  FIR  filter  time  needed  is  MT  which  is  also  the 
duration  of  its  transient  response.  l-'!'  .  j  define  this  time  to  be  also  the  settling 
time  of  a  recursive  filter  with  the  sa;i,e  specification  as  the  nonrecursive  filter. 

Even  though  longer  dwell  times  are  needed  for  a  filter  with  sharper  characteristics, 
it  is  not  the  filter  dwell  time  that  dictates  the  radar  update  rates.  These  are 
determined  by  the  desired  accuracy  of  moment  estimates.  For  instance,  with  a 
continuous  pulse  train  one  may  use  more  or  less  than  M  samples  (M  is  the  FIR  filter 
lenqth  or  length  of  an  equivalent  FIR  when  the  recursive  filter  is  used)  in  the 
autocovariance  processor,  and  the  velocity  or  spectrum  width  estimates  after 
filtering  may  be  updated  more  frequently  than  every  M  pulses.  Thus  one  can  always 
match  the  update  times  to  the  antenna  rotation  rate  (see  Appendix  D).  With  a  change 
in  rotation  rate  the  spectrum  width  of  ground  clutter  changes  and  to  maintain 
proper  filtering  the  notch  must  be  adjusted  according  to  (1.9).  For  instance  at 
6  rpm  the  passband  cut-off  velocity  should  be  v  =  2.5  -  3  m-s”1.  The  filter  with 
Vp  =  3  m*s"'  is  discussed  thoroughly  in  a  latersection  where  it  is  labeled  as 
filter  number  2. 
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2.  TRANSMITTED  WAVEFORMS,  SCANNING  STRATEGY  AND  CLUTTER  CANCELING 

Methods  for  canceling  clutter  are  intimately  tied  to  the  signal  design  and 
processing  schemes  of  the  weather  radar.  For  brevity,  we  consider  only  techniques 
that  are  economically  feasible  and  intuitively  satisfying.  Thus,  at  the  onset  we 
will  assume  that  an  8  or  more  bit  A/D  converter  is  available  so  fiat  its  saturation 
would  be  unlikely.  If  a  pulse  to  pulse  automatic  gain  control  circuit  (AGC)  is 
used,  it  should  have  a  large  time  constant  so  as  to  smooth  variations  presented  to 
the  filter.  Thus,  we  will  ignore  the  effects  of  AGC  since  these  in  principle  can 
be  either  accounted  for  or  avoided. 

For  precise  moment  estimation,  a  recursive  filter  (infinite  impulse 
response  HR)  or  a  finite  impulse  response  filter  (FIR)  are  feasible.  Delay  line 
cancel ers  cannot  be  employed  in  conjunction  with  the  pulse  pair  processor  (see 
Groginsky  and  Glover,  1980)  and  our  attempt  to  do  filtering  on  the  composite 
autocovariance  of  signal  plus  clutter  had  meager  success  (see  Appendix  C). 

We  consider  a  50  dB  rejection  in  the  stopband  and  a  1  dB  ripple  in  the  pass- 
band  reasonable  goals  to  achieve  since  a  recursive  filter  with  such  characteristics 
in  steady  state  has  already  been  designed  by  Raytheon  engineers  (Groginsky  and 
Glover,  1980).  Incidentally,  no  more  than  50  dB  of  cancellation  is  possible  if  an 
8-bit  A/D  converter  is  used  because  when  the  clutter  spans  the  full  8-bit  range, 
it's  quantization  noise  is  50  dB  below  its  average  power.  The  performance  of  this 
recursive  filter  on  a  contiguous  pulse  train  is  shown  in  the  next  section  while 
the  FIR  filter  and  the  number  of  pulses  needed  to  achieve  the  desired  characteris¬ 
tics  is  discussed  in  Appendix  B. 


£ 
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The  type  of  transmitted  signal  has  a  tremendous  bearing  on  which  of  the  two 
filters  would  yield  better  results.  To  take  full  advantage  of  the  recursive 
filter,  the  transmitted  sequence  must  be  uniform  and  uninterrupted.  This  would 
then  require  two  frequencies — one  for  velocity  measurements  and  the  other  for 
reflectivity.  The  performance  of  a  dual  frequency  radar  in  clutter  can  be 
exactly  predicted  from  the  knowledge  of  the  filter  type.  Also  the  required  update 
rates  to  achieve  tolerable  errors  in  estimates  can  be  obtained  from  published 
curves  (Zrnic' ,  1979).  We  will  not  dwell  on  this;  rather,  in  the  next  few  para¬ 
graphs,  we  discuss  an  interlaced  pulse  scheme  and  a  corresponding  scanning  strategy 
that  allows  considerable  clutter  cancelling  at  low  elevation  angles  where  it  counts. 

Let  us  assume  that  the  PRT.T  ,  of  the  radar  would  be  variable  and  somewhere 
between  the  limits  as  in  Table  2.1. 

Table  2. 1 


PRT  Ts 

Unambiguous  velocity  va 

a 

Unambiguous  range  r  for 
velocity  measurements 


.781  ms 

1.2  ms 

32  m-s"1 

20.8  m*s 

117  km 

180  km 

-1 


Ik  _i~ 


Operation  with  the  longest  PRT  produces  maximum  dwell  time,. and  if  we  do  not 
allow  more  than  50  percent  degradation  in  the  apparent  antenna  pattern  due  to 
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rotation  during  acquisition,  we  find  the  relationship  between  rotation  rate  a 
(°s~r),  T  and  the  number  of  samples  M  for  a  1°  antenna  beamwidth  (see  Appendix 
0)  s 

aMTs  =  1  (2.1) 


Both  the  FIR  filter  and  HR  require  about  64  samples  for  good  clutter 
canceling  (i.e.,  notch  width  2  vD  •  4  m-s-1 ,  and  40  to  50  dB  cancelation,  see 
Fig.  B.6).  Let's  accept  M=64  ana  let  us  add  another  8  PRT's  for  reflectivity 
estimation  as  in  Fig.  2.1. 


Figure  2.1  Interlaced  pulse  train  for  separate  velocity  and  reflectivity  estima¬ 
tion.  At  low  elevation  angles  about  64  pulses  are  needed  to  achieve  good  clutter 
canceling ,  while  at  higher  elevation  angles  M  could  he  16.  To  estimate  reflec¬ 
tivity  with  acceptable  accuracy 3  a  broadband  pulse  must  be  employed  in  conounc¬ 
tion  with  some  range  averaging  (Krehbiel  and  Brook ,  1979). 

The  pulse  for  reflectivity  would  need  to  be  a  broadband  signal  (perhaps  5  to 
10  MHz)  to  reduce  statistical  fluctuations  of  reflectivity  estimates  (Krehbiel  and 
Brook,  1979).  From  (2.1)  we  find  the  rotation  rate  a=12°s"  (or  2  rpm).  Such  a 
slow  rotation  rate  need  not  be  maintained  throughout  the  whole  scan.  Only  the  two 
lowest  elevation  angles  (.5°  and  1°)  would  definitely  require  that  much  canceling. 

At  higher  elevations,  the  cancel er  can  be  bypassed  and  the  rotation  rate  increased. 
For  instance,  16  contiguous  pulses  with  some  averaging  in  range  are  sufficient  to 
provide  adequate  accuracy  (Figs.  2,2  and  2.3).  With  8  more  PRT's  for  reflectivity  we 
find  M*=24  (where  M‘=is  the  number  of  system  periods,  T  ,  needed  for  both  reflec¬ 
tivity  and  velocity)  and  a*‘36°s_  or  6  rpm.  A  possible  scanning  strategy  employing 
the  above  principle  is  listed  in  Table  2.2. 


Table  2.2 


Elevation  angle 

M* 

Rotation  rate 

Time 

.5°,  1° 

72 

2  rpm 

1  min 

2°,  3°,..., 13° 

24 

6  rpm 

2  min 

15°,  17°,  19°,.. .,25° 

24 

6  rpm 

1  min 

While  this  may  not  be  an  optimum  strategy,  it  illustrates  that  4-minute 
update  times  are  realistic  with  significant  ground  clutter  canceling  capability  at 
the  two  lowest  elevation  angles. 

If  a  notch  width  2  v  of  8  m^s"^  is  acceptable,  then  a  32  point  weighted 
Fourier  transform  may  be  adequate  (see  Appendix  B),  in  which  case  the  rotation 
rate  at  the  lowest  elevation  angles  would  increase  by  a  factor  of  two.  This  would 
then  allow  canceling  at  two  additional  elevation  angles  (2°  and  3°).  Also  use  of 
initialization  on  a  recursive  filter  could  further  reduce  the  number  of  needed 
contiguous  pulses.  This  is  an  area  requiring  more  research. 
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SPECTRUM  WIDTH  <rv  (m  s'1) 
(T»  *  1.2  ms,  Vo  *  20.8  m  s*1) 


SPECTRlM  WIDTH  ct  y  (m  s'1) 
(Ts  =  78l^.s,  va  =  32  ms) 


Figure  2. 2  Standard  deviation  of  velocity  estimates  from  the  pulse  pair  algorithm 
with  16  pulses  (i.e.,  15  pairs).  The  upper  and  right  scales  are  for  Ts  *  1.2  ms; 
the  lower  and  left  scales  are  for  T  =  781  ys. 

o 

SPECTRUM  WIDTH  <rv  (m  s'1) 

(T,  =  1.2  ms,  v0  *  20.8  m  s'1) 


0  2  4  6  8  10  12  14 


SPECTRUM  WIDTH  <rv  (ms*») 
(T,=  781/is,  v0  *  32  m  $*') 


Figure  2.3  Standard  deviation  of  spectrum  width  estimates  from  the  pulse  pair 
algorithm  with  16  pulses. 


18 


1 

& 

s 


ii 


1 


3.  GROUND  CLUTTER  FILTER 

This  section  examines  practical  cancelers  that  could  function  in  conjunction 
with  the  pulse  pair  processor.  A  recursive  filter  that  operates  on  contiguous 
pulse  pairs  is  examined  first. 


3.1  Recursive  Filter-Contiguous  Pulse  Train 

The  recursive  filter  we  consider  is  a  three  pole  elliptic  filter  with  a 
50  dB  rejection  in  the  cutoff  band  and  a  transmission  zero  at  DC,  a  1  dB  ripple 
in  the  passband  and  a  variable  notch  width.  This  type  of  a  filter  has  been  built 
by  Raytheon  for  the  AFGL  (Groginsky  and  Glover,  1980).  In  Table  3.1  are  listed 
three  filters  that  we  have  designed  and  analyzed.  Figure  3.1  is  the  appropriate 
flow  chart  that  defines  the  K  coefficients. 


Table  3.1.  Coefficients  for  the  high  pass  elliptic  filter. 


Filter 

number 

fP 

VP 

fs 

vs 

K1 

k2 

K3 

k4 

Hz 

(m/s) 

Hz 

(m/s) 

1 

16 

0.8 

4.5 

0.225 

1.999621 

1.959927 

0.965686 

0.858737 

2 

60 

3.0 

17.0 

0.85 

1.994598 

1.809719 

0.895496 

0.455619 

3 

80 

4.0 

22.6 

1.13 

1.990548 

1.715192 

0.863963 

0.339933 

Figure  3.2.  Frequency  (Doppler  velocity)  response  characteristic  of  the 
recursive  Filter  1.  The  Nyquist  interval  is  65.1  m-s-1;  in  subsequent  simulations 
it  is  slightly  reduced  to  64  m-s~ 


For  illustration,  the  transfer  characteristic  of  the  Filter  1  is  shown  on 
Fig.  3.2,  which  demonstrates  that  the  specifications  have  been  met.  The  other 
two  filters  have  very  similar  characteristics  but  with  a  wider  notch. 

The  recursive  filter  operates  on  the  I,  Q  digital  video  samples  before  the 
pulse  pair  processor. (Fig.  3.3). 


Figure  3.3  Block  diagram  of  data  flow.  I  and  Q  are  the  digitized  video  signals 
the t.  must  be  filtered  prior  to  pulse  pair  processing . 


In  addition  to  the  mean  velocity  v  and  spectrum  width  °  estimates,  the  processor 
provides  an  estimate  of  echo  power  P.  v 

To  illustrate  the  performance  of  the  tandem  filter-pulse  pair  processor, 
simulations  on  synthetic  time  series  data  and  on  some  real  data  were  conducted. 
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Roundoff  effects  on  coefficients  and  finite  word  lengths  of  data  were  not  considered. 
Groginsky  and  G’over  (1980)  show  that  coefficients  with  12  to  16-bit  word  lengths 
(depending  on  radar  PRT)  do  not  degrade  the  filter  performance.  Also,  data  with 
8  or  more  bits  would  have  ground  clutter  quantization  noise  more  than  50  dB  below 
the  ground  clutter  power  when  it  spans  the  whole  A/D  range. 

On  Fig.  3.4a,  b  are  the  spectra  obtained  from  real  time  series  data  before 
and  after  the  recursive  Filter  #2.  To  reduce  transient  effects,  the  first  40 
points  of  128  time  series  data  were  tapered  with  the  raised  cosine.  The  power 
spectrum  of  the  last  64  points  is  taken.  In  either  case  prior  to  taking  the 
spectrum,  data  were  multiplied  with  the  von  Hann  window.  It  is  evident  from  the 
comparison  of  3.4a  and  b  that  a  50  dB  notch  was  achieved  and  that  the  weather  peak 
(at  -7  m-s-  )  that  was  15  dB  below  clutter  would,  after  filtering,  produce  a 
viable  mean  velocity  and  width  estimates.  Note  also  that  if  the  DC  line  and  one 
spectral  power  either  side  of  it  are  removed  from  the  von  Hann  weighed  periooogram 
(Fig.  3.4a),  the  weather  echo  would  stick  12.5  dB  above  the  clutter  peak.  After 
the  recursive  filter,  the  weather  peak  is  22  dB  above  the  highest  clutter  residual 
which  means  that  at  least  37  dB  of  clutter  power  has  been  canceled.  This  filter's 
stopband  cutoff  is  between  -.85  m-s“  and  .85  m-s-  where  we  see  over  50  d^  of 
cancellation  was  achieved  (Fig.  3.4b).  Passband  cutoff  is  between  -3  m-s"  and 
3  m-s-  and  the  residual's  peak  is  at  4  m-s-  where  the  slight  overshoot  of  this 
filter’s  frequency  response  occurs. 

Digital  computer  simulation  cf  the  block  diagram  on  Fig.  3.3  was  done  to  test 
the  quality  of  the  output  after  filtering.  The  input  time  series  data  were 
synthesized  such  that  they  contained  a  clutter  and  weather  signal  of  Gaussian 
spectral  shapes  and  had  a  Gaussian  statistics  (see  Zrnic',  1975,  for  details  and  -| 
Appendix  E  for  program  flow  charts).  Unambiguous  velocity,  v  ,  was  set  to  32  m-s-  . 
In  the  results  that  follow,  the  filter  was  allowed  to  reach  steady  statej'.efore  64 
points  were  processed  with  the  pulse  pair  estimator  to  obtain  P,  v,  and  £- 
Seventy  of  the  estimates  were  averaged  to  obtain  the  mean  curves  that  arevshown  on 
Figs.  3.5  through  3.16.  Graphs  on  Fig.  3.5  suggpst  that  filters  with  wider  notch 
(numbers  2  and  3)  overall  create  less  bias  in’ the  power  estimate  P  than  the 
narrowest  filter,  except  near  the  zeroth  velocity.  Also  evident  is  the  need  for  a 
very  good  match  betw.Qn  the  notch  width  and  the  width  of  the  clutter  spectrum.  It 
seems  that  filter  1  would  be  better  matched  to  a  clutter  width  of  .1  m-s  rather 
than  the  .2  m-s"  used  as  example  here.  We  point  out  that  estimation  of  power 
from  the  coherent  data  (I,  Q)  would  be  attempted  only  if  no ‘echoes  are  overlaid. 
Otherwise,  the  power  should  be  estimated  from  the  reflectivity  channel.  Modest 
canceling  of  ground  clutter  in  the  reflectivity  channel  is  possible  (see  Preprints 
of  the  19th  Radar  Meteorology  Conference)  but  will  not  be  discussed  here. 

Graphs  of  the  estimated  meai  velocity  versus  true  mean  velocity  for  the  3 
filters  are  shown  in  Figs.  3.6,  3.7,  and  3.8.  These  (and  other  figures  not 
shown)  suggest  that  the  formula  (1.9)  for  the  notch  width  is  essentially  correct. 
This  formula  states  that  the  cauceler  will  perform  well  if  clutter  widths  are  less 
than  v  /4,8.  For  instance.  Filter  1  is  matched  to  the  clutter  width  of  = 

.16  m-s  "  ;  for  this  width  (see  figures)  and  smaller  widths  (not  shown)  ivperforms 
wel4  but  we  see  that  at  a  =  .2  m-s"  the  estimate  is  good  only  up  to  20  dB  of  C/S 
ratio.  Filter  2  is  matched  to  a,  =  .63  m-s"  and  we  note  that  the  performance  is 
good  for  a  =  .2  and  .5  m-s"  and  is  very, degraded  at  a  =  1  m-s"  as  it  should 
be.  Filter  3  is  matched  to  o  =  .83  m-s  ,  and  we  notecbeginning  of  degradation 
at  a  =  1  m-s.  Anoverbiasid  estimate  (by  2  m-s"  or  less)  occurs  at  large 
signal  widths  (8  m-s"  )  and  good  C/S  ratios  (Fig.  3.8)  due  to  canceling  of  low 
velocity  components  in  the  weather  spectrum. 
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Figure  3.4.  a)  Power  spectrum  of  time  series  containing  ground  clutter  and 
weather  signal.  Logarithms  of  powers  with  respect  to  the  peak  are  plotted 
(X's)  and  the  solid  curve  is  a  five-point  running  average  for  visual  clarity, 
b)  Same  as  a)  but  time  series  data  were  filtered  with  the  elliptic  recursive 
Filter  #2. 


■t 


•a 


22 


t  it  't*  : 
.{**0  & .16 
'Ar+tl  .0 


nut*  *. 

Cl*C  =  .2 

S!K£  :1  *0 


4  S  \2  15  2D  24  20  32 

TRUE  ftftN  VELOCITY 


FILTER  1 
COO  £.S 
5100  si  0 


hi'kl 
too  i.tl 
5100  5.3 


■ft—T'V  \Y“iT  a  ‘ 

tRi>  *tfcGN  VM  L:.:fY 


clLT  *  1 
CUO  3.2 
5100  =C-0 


%  4  8  12  16  2C  24  20  32 

TR'€  rtEflN  VELOCITY 


FILTER  1 
CLIO  s.$ 
SlOO  *8.0 


8  12  16  20  24* 

TRUE  rrffiN  vaxiTv 


filter  i 

COO  £l-C 

s:oo  si.o 


8  12  16  20  24* ‘*8  32 

TRUE  «E«N  VELOCITY 


filter  i 

Cue  :)  .0 
SlOO  *8  0 


6  12  16  20 
TRUE  nEttf  VELOCITY 


0  12  16  20 
TRUE  It  AN  VELOCITY 


Figure  3.  6.  Estimated  mean  velocity  versus  true  mean  velocity.  Signal  widths  are 
1  m*s~l  (on  the  left)  and  8  m*  s~l  on  the  right.  Clutter  width  increases  from 
top  to  bottom j  and  is  indicated  togetiier  with  clutter  to  signal  ratio. 


The  filters  perform  similarly  for  the  whole  range  of  weather  signal  widths 
between  1  and  8  m-s"  .  To  further  emphasize  the  point,  larger  C/S  ratios  are 
shown  on  Fig.  3.9.  The  sensitivity  to  clutter  width  is  apparent.  Note  that  40  dB 
of  C/S  can  be  tolerated  for  widths  smaller  than  the  matching  width  (0.83  m-s  ).■ 
Although  the  cancel er  removes  50  dB  of  clutter  in  the  spectrum  domain,  the  pulse 
pair  processor  operates  well  when  spurious  peaks  are  about  10  dB  below  the  signal 
peak,  and  therefore  50  dB  of  C/S  ratio  (per  pulse)  results  in  biased  estimate 
(middle  of  Fig.  3.9).  The  bias  is  quite  small  for  oc  =  0.2  m-s"',  but  that  is  be¬ 
cause  at  such  narrow  clutter  spectrum  width  more  than  50  dB  cancelation  is  achieved 
by  the  steep  region  of  the  filter  around  zero  (see  Fig.  3.2  for  the  filter  character¬ 
istic). 

Behavior  of  the  estimates  around  zero  velocity  is  of  interest  since  that  is 
where  the  filter  characteristic  changes  most.  For  the  wide  notch  Filter  #3 
(2v  =8  m-s"  ),  the  mean  velocities  are  plotted  on  Fig.  3.10.  The  estimates  aife 
rea§onably  well -.behaved  except  for  a  rather  anomalous  curve  when  a  -  .2  m-s" 
and  a  =  1  m-s"  .  It  is  not  clear  if  this  anomaly  is  real  or  computational. 

Apartvfrom  it,  the  estimates  are  biased,  but  this  should  not  present  a  problem  for 
meteorological  interpretations  since  biases  are  small  and  predictable.  With  the 
narrower  filter,  the  transition  through  zero  is  much  smoother  (Fig.  3.11). 

From  the  above  consideration,  it  is  essential  to  have  the  passband  cutoff  v 
equal  to  or  slightly  larger  than  4.8  a  .  The  filter  was  operated  in  steady  state 
and  provided  64  samples  to  the  pulse  pair  estimator.  Under  these  conditions, 
excellent  mean  velocity  estimates  are  obtained. 

Conclusions  concerning  spectrum  width  estimates  (Figs.  3.12,  3.13,  3.14,  and 
3.15)  are  very  similar  to  the  mean  velocity  estimates.  However,  the  bias  can  be 
worse  under  some  circumstances.  For  instance.  Filter  2  is  matched  to  clutter 
width  of  0.63  m-s"  ; thus  at  smaller  widths  it  should  perform  well.  In_l[ig.  3.13 
we  see  that  at  a  mean  velocity  of  24  m.s"  and  clutter  width  of  .5  m-s  the 
signal  width  estimate  has  a  significant  bias  even  though  the  clutter  to  signal 
ratio  is  only  30  dB.  Behavior  of  width  estimates  in  the  transition  region  is 
plotted  on  Figs.  3.16a  and  b.  Biases  at  mean  velocities  below  2  m-s-1  are  evident 
when  data  are  filtered  with  the  wide  notch  Filter  3,  but  such  cases  can  be  auto¬ 
matically  recognized. 

3.2  Recursive  Filter,  Stepped  Input  and  Initialization 

The  main  disadvantage  of  the  recursive  elliptic  filter  is  its  transient 
response.  Although  the  steady  state  frequency  response  of  HR  filters  is  similar 
to  comparable  FIR  filters,  infinite  duration  of  their  impulse  response  has  limited 
their  use  in  radar  application.  At  low  elevation  angles  (<  0.5°)  the  clutter  power  is 
several  orders  of  magnitude  stronger  than  the  signal  power.  Thus,  when  interlaced 
samples  (batch)  for  velocity  and  reflectivity  estimation  are  used,  the  clutter 
appears  as  a  step  to  the  filter  which  produces  severe  ringing  at  the  output  and  in 
effect  overwhelms  the  signal.  To  reduce  this  effect,  either  enough  time  must  be 
allowed  for  the  transients  to  settle  before  pulse  pair  processing  is  attempted  or 
some  other  means  of  their  suppression  must  be  employed.  Because  settling  time  is 
tens  of  pulses  (Fig.  3.17),  we  consider  the  alternative  which  is  initialization 
(Fletcher  and  Burlage,  1972).  The  purpose  of  initialization  is  to  load 
the  memory  elements  of  the  recursive  filter  with  their  anticipated  steady 
state  values.  These  values  can  be  estimated  from  the  first  incoming 


TRUE  MEAN  VELOCITY 


MUM  » 
awo  *.>6 
VtUXtl*  tH.0 


f  HIM  i 
a*c  «.is 
vtlDCHT  *120 


rtUM  i 
cue  *. $ 
veto cm  *2i.o 


2  4 

TRj£  HERN  WIO!* 


4  6 

true  ni«N  width 


2  4  6 

TRl 1  OtRN  WIDTH 


HUM  • 

C.IC  '  2 

*tiocr '  *•«  o 


MUM  \ 
a  >o «  z 
vtuxm  »?*  o^. 


'*jl  hern  *:3»h 


IRlt  HCflN  WIDTH 


TRUE  HERN  H10TH 


HUM  J 
a»c  *  s 

vtaxm  »*  <3 


•k*  hUiTm 


IRV*  PfWS  Wl.MW 


2  4  I 

TRUt  ntftN  WIDTH 


♦  tirw  : 
aw)  ii.o 

VUOCJIT  t?»  0 


TRUE  ntflN  WIDTH 


Figure  S.12.  Estimated  mean  signal  spectrum  width  versus  true  width  for  simulated 
time  series  data  processed  by  the  pulse  pair  processor.  Filter  1  is  used  and 
clutter  width  increases  from  top  to  bottom.  The  mean  velocity  of  signal  spectrum 
is  4  m* s~^  (left  graphs),  12  m'S~^  (inside  graphs ),  and  24  m' s~^  (right  graphs). 
Ground  clutter  to  signal  power  ratios  are  indicated  and  range  from  0  to  SO  dB. 
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3.13.  Scone  as  3.12  but  filter  2  is  used  and  clutter  width  of  .16  m‘S~ ^  is  not  ?| 
idered.  ! 
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Figure  3.16a.  Estimated  moan  width  versus  true  mean  width  when  the  true  velocity  J 

is  in  the  notch  of  filter  3.  The  mean  velocity  increases  in  steps  of  1  m-  s~7  ;; 

from  left  to  right.  1 


Figure  3.16b.  Same  as  a)  but  filter  1  was  used ,  and  two  mean  velocities  .1  and 
3.5  m*  s~*  are  considered.  Other  velocities  between  these  two  produce  no 
noticeable  difference  in  the  results. 
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pulse.  When  the  ground  clutter  is  a  perfect  D.C.  line  and  no  weather  signal  is 
present,  the  initialization  will  bring  immediately  the  filter  into  steady  state. 

Alas,  the  clutter  has  finite  width  and  superposed  to  it  is  the  weather  signal 
which  together  degrade  the  initialization. 

_1  For  the  third  order  filter  on  Fig.  3.1,  the  output  of  the  first  memory  element, 

Z  ,  would  be  set  to  A,/(i-Kn)  where  A,  is  the  digital  value  of  the  first  sample. 

The  outputs  of  the  othdr  twoVlements  must  be  set  to  0.  With  this  the  output 
after  the  first  time  sample.  A,,  is  0  and  after  the  second  one,  A2,  it  is  A^-A,. 

Thus  the  filter  immediately  acts  as  a  delay  line  canceler  so  that^a  large  D.C. 
value  does  not  get  through.  Simulation  of  initialization  on  a  digital  computer 
is  briefly  described  in  Appendix  E,  subroutine  ELLIPSE  1  • 

Quantitative  measure  of  the  performance  of  initialization  was  obtained  from 
simulations.  First,  we  will  discuss  transient  effects  on  the  moment  estimates 
with  no  initialization.  Figure  3.18a  shows  plots  of  the  estimated  mean  velocity 
versus  true  mean  velocity  for  the  case  of  no  clutter  and  C/S  =  10  and  20  dB  with¬ 
out  any  initialization  or  smoothing.  Blocks  of  64  time  samples  were  applied  repeatedly 
(70  times)  to  the  recursive  filter,  and  the  estimates  of  the  mean  velocity  and  spectrum 
width  were  obtained  from  a  64  sample  pulse  pair  algorithm.  The  estimates  we re  then 
averaged  to  generate  the  mean  yalues.  Throughout  this  simulation,  the  clutter  spectrum 
width  was  set  to  a  =0-16  m-s"  which  is  the  upper  limit  for  cancellation  with 
Filter  No.  1.  An  unbiased  estimate  of  the  mean  velocity  is  obtained  in  the  absence 
of  flutter.  Bias  is  quite  evident  for  both  10  and  20  dB  C/S  ratios  and  typical 
signal  spectrum  widths  1.0  m-s"  <  a  <8.0  m-s  .  The  width  estimate  has  very 
little  bias  when  there  is  no  clutter  present  (Fig.  3.18b).  For  C/S  =  10  dB  and 
velocities  12.0  m-s,  the  width  estimate  is  very  poor  and  is  further  degraded 
at  C/S  =  20  dB. 

These  biases  are  mainly  there  because  the  filter  was  in  transient  during 
approximately  the  time  half  of  the  samples  were  applied.  Obviously  these  samples 
were  not  properly  filtered  and  their  residuals  have  created  biases  in  the  mean 
velocity  and  spectrum  width. 

To  suppress  the  transient  reponse,  initialization  technique  was  used,  i.e., 
the  first  three  memory  elements  were  set  to  anticipated  steady  state  values. 

Figure  3.19a  shows  the  plot  of  estimated  mean  velocity  versus  true  mean  velocity 
for  C/S  10  and  20  dB.  The  estimated  values  have  very  little  bias  for  signal 
spectrum  widths  ov*1.0  m-s.  An,unbiased  estimate  of  the  width  (Fig.  3.19b)  is 
obtained  for  velocities  8.0  m-s"  and  C/S  =  10  dB.  At  C/S  =  20  dB  the  bias 
associated  with  the, estimated  signal  spectrum  width  is  large  except  when  the  mean 
velocity  is  4  m-s"  .  In  fact,  the  results  are  better,  i.e.,  biases  are  one-half 
of  the  ones  when  C/S=10  dB  and  no  initialization. is  applied  (Figs.  3.18a  and  b).  Thus, 
in  conclusion,  filtering  over  64  pulses  with  no  initialization  is  equivalent  to 
about  15  dB  S/C  improvement,  and  initialization  adds  another  15  dB  of  improvement. 

In  other  words,  without  any  filtering  a  10  dB  S/C  ratio  for  velocity  and  15  dB  for 
spectrum  width  is  needed  to  accurately  estimate  moments.  With  recursive  filtering 
of  stepped  inputs,  64  pulses  long  at  C/S  of  0  -  5  dB  one  may  obtain  valid  v  and 
a  while  initialization  further  raises  this  ratio  to  15-20  dB.  In  conclusion, 
recursive  filtering  on  interlaced  samples  implies  a  15-20  dB  penalty  in  ground 
clutter  rejection  as  opposed  to  filtering  a  contiguous  pulse  train  in  steady  state. 
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Figure  3.18.  a)  Mean  velocity  estimates  frm  the  pulse  pair  algorithm  when  a 
64  sample  pulse  train  is  filtered  by  a  recursive  filter  No.  1  and  no  provision 
to  suppress  transients  is  applied,  b)  Mean  spectrum  width  estimates  under 
same  conditions  as  in  a). 
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It  may  be  possible  to  reduce  this  penalty  by  ignoring  the  first  few  samples 
at  the  filter's  output  (after  initialization  and  processing  the  remaining  ones. 

At  present  it  is  not  clear  what  is  the  minimum  number  of  samples  that  would  allow 
good  clutter  canceling  and  how  many  of  the  first  samples  (if  any  after  initializa¬ 
tion)  should  not  be  processed. 

We  have  determined  that  initialization  is  not  detrimental  to  the  weather 
signal  even  if  ground  clutter  is  weak.  For  instance,  a  64  point_ spectrum  of  real 
signal  and  clutter  is  shown  in  Fig.  3.20a.  The  signal  at  -8  m*s"  is  10  dB  above 
clutter.  Next  to  this  graph  (Fig.  3.20b)  is  the  spectrum  of  the  filtered  signal 
without  initialization,  i.e.,  the  64  time  samples  were  passed  through  the  eliptic 
filter.  On  Fig.  3.20c  is  the  spectrum  of  the  same  time  samples,  but  the  recursive 
filter  has  been  initialized  and  on  3.20d  the  same  time  samples  have  been  used 
except  that  previous  64  samples  were  also  passed  through  the  filter  to  make  it 
reach  steady  state.  We  note  that  the  clutter  has  been  suppressed  to  about  45  dB 
below  signal  when  the  stepped  input  is  applied  to  the  filter.  The  residual  is 
50-55  dB  below  signal  when  initialization  is  applied  and  is  even  further  down  when 
the  filter  has  reached  steady  state.  The  weather  spectrum  is  not  affected  by  the 
three  variations  in  the  filtering. 

Very  similar  conclusions  are  reached  if  the  clutter  power  is  larger  than  the 
signal  power.  An  example  on  Fig.  3.21a  shows  a  clutter  to  signal  ratio  of  about 
12  dB.  Filtering  with  no  initialization  removes  considerable  clutter  (Fig.  3.21b), 
but  initialization  further  decreases  the  residual  by  5  to  10  dB  (Fig.  3.21c). 

Of  course,  the  smallest  residual  is  obtained  if  the  filter  reaches  steady  state, 
which  it  approximately  does  when  one  also  applies  the  preceding  64  time  samples 
(Fig.  3 . 21 d ) .  In  this  example,  a  fairly  large  image  (about  12  dB  below  the  peak 
of  the  weather  signal  is  most  likely  due  to  amplitude  imbalance  of  the  inphase 
and  quadrature  components. 


3.3  Other  Techniques 

Two  more  approaches  for  filtering  ground  clutter  out  of  an  interlaced  pulse 
sequence  were  considered.  In  one,  time  varying  recursive  filters  with  periodic 
coefficients  meant  to  compensate  for  the  unequal  pulse  spacing  were  tried.  However, 
the  introduced  periodicities  created  strong  spectral  peaks  (only  5  dB  below  the  DC 
level)  which  could  not  be  eliminated.  A  finite  impulse  response  filter  with  time 
varying  coefficients  would  give  a  better  frequency  response  but  due  to  large 
storage  requirement  and  the  need  to  change  weights  at  every  step  makes  this 
filter  impractical . 

Still  another  possibility  is  to  use  the  recursive  filter  at  a  reduced  rate- 
let's  say  every  kT  where  k  is  a  small  integer.  But  this  introduces  transmission 
zeroes  at  frequencies  (n=l ,2,, . . ,k)  which  would  eliminate  considerable  amount  of 
weather  signal.  Furthermore,  filtering  of  pulses  between  the  first  and  the  k'th 
is  not  trivial . 
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Figure  3.21.  a)  Spectrum  Of  weather  and  clutter  for  C/S  »  12  dB. 
b)  Time  series  that  generated  a)  was  filtered. 

a)  Time  series  was  filtered  and  the  filter  was  properly  initialized. 

d)  128  time  samples  were  passed  through  the  filter  and  the  last  64  corresponding 

to  the  ones  in  a)}  b)3  and  c)  are  Fourier  transformed. 
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4.  CONCLUSIONS 

Study  of  the  ground  clutter  problem  at  the  Norman  site  indicates  that  clutter 
cross  sections  of  40  to  50  dB  below  1  nr/nr  are  typical.  These  values  are  small 
and  the  clutter  extends  in  range  to  about  30  km.  If  a  rejection  in  a  stopband  of  50  dB 
can  be  achieved,  then  gust  fronts  with  at  least  -10  dBZ  reflectivity  would  be  observed 
in  all  but  400  km2  area  at  0°  in  elevation.  Thus,  only  in  25  percent  of  the  clutter 
area  such  weak  signals  could  not  be  retrieved.  On  the  other  hand,  reflectivities 
of  20-30  dBZ  would  produce  viable  spectral  moments  everywhere  except  at  a  few 
spots.  Thus,  a  similar  site  would  be  clutter  free  for  storm  observations. 

Doppler  spectrum  width,  o  of  the  clutter  signal  was  between  &1  and  0.5  m-s~T 
with  the  mean  value  of  .25  m.sc_l  at  an  antenna  rotation  rate  of  10°  s~  .  The 
design  of  clutter  filters  is  very  much  dependent  on  this  width  and  for  proper 
cancellation  the  passband  » utoff  velocity  Vp  must  be  larger  or  equal  to  4.8  ac. 

Several  canceling  schemes  were  investigated,  and  the  most  efficient  one  is  with 
a  third  order  recursive  filter.  This  filter  achieves  a  50  dB  rejection  in  the 
stop  band  with  total  annihilation  of  DC.  One  dB  ripple  is  in  the  passband,  and 
the  ratio  of  passband  cutoff  v  to  stop  band  cutoff  v  is  about  3.5.  The  filter 
operates  best  in  steady  state, pbut  it  can  also  be  madi  to  operate  in  transient  by 
properly  initializing  its  memory  elements.  Performance  on  64  simulated  and  real 
time  samples  show  that  no  more  than  10  dB  loss  in  effective  canceling  over  steady 
state  occurs  if  the  filter  is  initialized.  Thus,  there  is  no  doubt  that  canceling 
will  be  effective  on  a  contiguous  64  pulse  sequence.  The  remaining  question  is, 
how  well  would  the  scheme  work  with  a  smaller  number  of  pulses  (16  or  32)  and 
if  the  first  few  pulses  from  the  filter's  output  should  not  be  processed.  In  our 
initialization  scheme  only  the  first  incoming  pulse  was  used  to  set  the  filter 
memory  elements,  after  which  the  system  was  left  to  itself.  One  may  perhaps  bene¬ 
fit  if  some  of  the  storage  elements  are  reinitialized  after  the  second  and  also 
after  the  third  pulse. 

The  autocovariance  processor  of  mean  velocity  requires  10  dB  of  signal 
to  clutter  power  ratio  while  for  second  moment  estimation  15  dB  insures  adequate 
accuracy.  Therefore,  the  two  estimates  will  not  behave  in  a  same  manner  for  a 
fixed  amount  of  clutter  canceling.  For  instance,  a  properly  initialized  filter 
would  allow  processing  of  samples  with  20  to  40  dB  of  C/S  for  mean  velocity  and 
15  to  35  dB  for  Doppler  spectrum  width. 

It  is  shown  that  a  scan  sequence  can  be  accomplished  in  about  4  minutes  if  the 
rotation  at  the  two  lowest  elevations  is  slow  (2  rpm)  and  at  the  next  18  levels 
it  is  increased  to  6  rpm.  The  purpose  here  was  not  to  investigate  scanning 
strategies,  but  rather  to  point  out  that  4  to  5  minutes  update  times  are  feasible 
with  good  clutter  canceling  at  low  levels  where  it  counts.  Other  possible  schemes 
such  as  interlaced  elevation  angles,  etc.,  would  also  need  to  have  antenna  speeds 
reduced  to  2  rpm  at  the  lowest  elevations.  Beside  the  variable  antenna  speed  it 
may  be  wise  to  consider  an  adaptive  canceling  scheme.  For  instance,  canceling 
should  be  done  only  where  needed,  i.e.  up  to  a  maximum  range  that  depends  on 
elevation.  Also  the  filter  notch  may  change  wUh  elevation. 


We  stress  that  clutter  canceling  throughout  the  volume  would  be  easiest  if  a  | 

dual  frequency  transmitter  is  employed,  because  then  the  recursive  filter  would  not  | 

be  subjected  to  transients  (if  block  AGC  is  utilized,  transients  due  to  gain  change 
are  easier  to  control).  Furthermore,  rotation  rate  at  the  lowest  elevations  could 
be  somewhat  Increased,  and  there  would  be  more  samples  for  velocity  estimation  at  | 


higher  elevations.  Also,  more  samples  for  reflectivity  calculations  would  be  made 
available.  This  last  parameter  may  prove  difficult  to  estimate  in  an  interlaced 
mode  with  a  long  string  of  velocity  pulses.  It  may  even  be  necessary  to  transmit  a 
broadband  pulse  for  reflectivity  estimation. 

The  part  of  weather  signal  that  is  within  the  stop  band  cutoff  of  the  filter 
is  totally  lost  but  that  is  a  small  penalty  since  zero  radial  velocities 
occupy  a  small  portion  of  a  storm  and  can  usually  be  identified.  But  large  spectrum 
widths  indicative  of  turbulence  and  other  hazards  would  still  be  measurable. 
Reflectivities  would  be  underestimated  if  the  signal  is  directly  in  the  notch,  and 
would  be  in  error  by  about  1  dB  in  the  immediate  vicinity  of  the  notch.  Therefore, 
canceling  should  be  done  selectively  at  those  ranges  and  elevation  angles  where 
ground  clutter  is  significant. 

Filtering  via  a  weighted  discrete  Fourier  transform  was  briefly  investigated 
(Appendix  B)  and  although  the  number  of  pulses  needed  to  achieve  comparable  charac¬ 
teristics  is  about  70,  the  filter  itself  is  not  subject  to  transients.  Therefore, 
its  effects  can  be  exactly  predicted  and  with  some  fitting  of  weather  and  clutter 
spectrum  models  it  may  be  even  possible  to  retrieve  considerable  information  about 
the  weather  signal  of  zero  velocity.  Another  possible  advantage  of  the  Fourier 
method  is  that  it  can  be  implemented  in  conjunction  with  a  phase  diversity  radar. 

By  this  we  mean  that  the  Fourier  transform  may  be  performed  on  samples  that  are 
phase  adjusted  for  the  first  trip.  Then  after  rejecting  clutter  spectrum  coeffi¬ 
cients,  an  inverse  transform  is  made  and  the  resulting  time  samples  can  be  phase 
readjusted  for  the  second  trip.  Another  transform  or  autocovariance  processing  of 
these  samples  will  reveal  the  echoes  from  the  second  trip  while  the  clutter  residual 
would  be  spread  uniformly  through  the  Nyquist  interval.  All  this  is  computationally 
very  intense;  It  is  not  clear  if  and  how  an  analogous  algorithm  with  the  recursive 
filter  could  unscramble  the  second  trip  echoes  in  a  phase  diversity  radar. 

The  first  trip  clutter  overlaid  on  storms  of  interest  in  the  second  trip  is 
very  bothersome  because  the  obscuration  area  increases  10  or  more  times  and  the 
second  trip  signal  strength  has  a  10  to  15  dB  disadvantage  due  to  range.  Effective 
ground  clutter  canceling  is  therefore  very  important  in  operational  Doppler  weather 
radar. 
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APPENDIX  A 


Errors  in  Pulse  Pair  Derived  Velocity  and 
Spectrum  Width  Due  to  Anomalous  Signal 


The  signal  of  interest,  with  power  'S'  and  the  anomalous  signal  (for  instance 
ground  clutter)  with  power  *C *  are  both  assumed  to  have  Gaussian  shaped  power 
spectra  with  mean  frequencies  f  and  f  and  widths  o  and  a .  The  error  we  are 
interested  in  is  the  difference5between  v  and  v  where  v^  is  the  pulse  pair 

velocity  for  the  total  signal  (signal  +  cTutter)  and  v  is  thi  pulse  pair  velocity 
of  the  signal  of  interest.  5 
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It  can  be  shown  that  e  is  a  function  of  |v  -  v  |,  i.e.,  the  absolute  value  of 
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Figure  A.l  shows  e  for  various  power  ratios.  For  typical  clutter  width  of 
0.1  m.s  to  0.5  m's"  and  signal  width  of  1  to  8  m>s  ,  the  ratio  p  /p  is  less 
than  1.3  dB  (Fig.  A. 2)  which  means  that  the  performance  of  the  pulsespaTr  mean 
velocity  (through  3)  depends  primarily  on  the  S/C  ratio.  From  Fig.  A.l  we  note 
that  a  signal  10  dB  larger  from  clutter  will  suffer  less  than  1  m-s’  error  in 
its  velocity  estimate.  Thus,  with  no  other  suppression,  a  10  dB  margin  of  signal 
over  clutter  can  be  considered  more  than  adequate  for  precise  Doppler  estimation 
via  a  pulse  pair  algorithm. 

Spectrum  width  error 

Formulas  for  the  spectrum  width  when  two  narrow  band  Gaussian  signal 
are  superimposed  can  be  found  in  Zrnic*  (1979).  The  formula  reads: 
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where  o.  is  the  biased  width,  v  the  mean  velocity  of  weather  signal  and  the  other 
symbols  are  defined  in  the  preceding  section.  The  plot  of  this  equation 
(Fig.  A. 3  and  A. 4)  shows  large  biases  for  narrow  widths,  which  are  also  strongly 
dependent  on  the  signal's  mean  velocity.  Doppler  spectrum  width  of  4  m*s”  indicates 
moderate  turbulence,  thus  if  we  want  to  be  near  or  below  that  value  for  typical 
widths  of  1-3  m*s  ,  we  need  a  margin  of  at  least  15  dB  between  signal  and  clutter. 
Results  of  simulations  on  the  same  figures  and  also  in  Appendix  C  confirm  the 
validity  of  these  curves. 
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Figure  A.l.  Magnitude  of  the  normalised  velocity  error  versus  the  magnitude 
of  the  difference  between  the  signal  mean  velocity  and  ground  clutter  mean 
velocity.  For  ground  based  radars  v  =  0.  Signal  to  clutter  power  ratio  is 
indicated >  and  the  curves  are  strictly  correct  for  a  <  av  «  v . 
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Figure  A.  2.  The  ratio  of  signal  and 
clutter  correlation  coefficients  for 
Gaussian  shaped  spectra.  Clutter 
widths  are  indicated.  Unambiguous 
velocity  v.  *  32  m-s~ *. 
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figure  A. 3.  Biased  spectrum  width  when 
signal  and  clutter  are  present.  The 
clutter  spectrum  width  is  less  than 
.5  m- s~l  and  unambiguous  velocity  is 
32  m‘S~l.  Mean  velocities  are  indi¬ 
cated  and  results  cf  simulations  for 
v=24  m-s~l  are  superimposed.  Dashed 
line  is  at  45°.  To  obtain  values  for 
any  other  unambiguous  velocity  va, 
divide  the  scale  by  32  and  multiply 
with  va. 
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APPENDIX  B 

Finite  Impulse  Response  -  FIR  Filter 


The  objective  is  to  estimate  the  number  of  sample  points  for  a  ground  clutter 
filter  with  a  1  dB  ripple  in  the  passband  and  a  stopband  rejection  of  -50  dB  (with 
respect  to  passband)  for  a  desired  set  of  passband  and  stopband  cutoff  frequencies. 
Figure  B.la  shows  the  amplitude  response  of  a  typical  FIR  high  pass  filter  with 
stopband  cutoff  radian  frequency  uis;  passband  cutoff  frequency  u>D;  passband 
deviation  $y,  and  stopband  deviation  52.  Figure  B.lb  shows  the  Magnitude 
response  described  in  terms  of  passband  parameter  e  and  stopband  parameter  A. 
Comparing  Figs.  B.la  and  B.lb,  e  and  A  are  related  to  5-j  and  62  as 


Figure  B.l.  Terminology  used  to  describe  high  pass  FIR  characteristics, 
a)  Amplitude  response ,  and  b)  Magnitude  response. 
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To  obtain  the  estimate  of  the  value  of  M  required  to  meet  specifications  on 
6-j,  <S2»  fp.  and  f$,  the  following  equation  can  be  used  (Rabiner  et  al.  1974). 
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f(6r  S2)  =  11 .01217+0-61244  log^  -U. 51244  log^ 
and  the  normalized  transition  width  Af  is  defined  as 
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Table  B.l.  The  number  of  sample  points  needed  for  various  passband  and 
stopband  cutoff  frequencies ,  j-4  1302  Hz  fva  =  32.5  m-s"1). 
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Figure  B.2  shows  the  number  of  sample  points  versus  the  notch  width  associated 
with  the  FIR  filter  where  the  width  is  defined  as  twice  the  passband  cutoff 
velocity.  From  the  figure  it  is  clear  that  in  order  to  achieve  narrow  widths,  one 
needs  to  collect  a  lot  of  sample  points  (M >125) ,  whereas  for  wider  widths  fewer 
samples  are  required  (M<75).  Therefore,  for  the  actual  design  of  the  filter, 
specifications  as  well  as  the  limitations  on  the  data  collection  must  be  considered. 
Because  the  pulse  repetition  time  for  these  calculations  is  768  ys,  the  dwell  time 
to  achieve  a  desired  cutoff  width  can  be  calculated  by  multiplying  M  on  Fig.  B.2 
with  768  ys.  This  is  shown  on  the  auxiliary  axis.  It  can  be  seen  by  evaluating 
(B.6)  for  this  type  of  filter  that: 


M  =  Dot  (6,  ,62}/Af 


or  from  (B.7c) 


MTS  *  Dro  ( 6-j ,  «2 )/  (fp  -  f$) 


(B.8) 


Thus  it  takes  a  minimum  dwell  time  MTs  to  achieve  desired  specifications. 

Two  hardware  implementations  of  this  filter  are  possible.  In  one  filter 
weights  are  determined  and  a  weighting  average  with  these  weights  is  computed,  i.e.,: 

k 

y(k)  =  l  w(n)x(k-n)  (B.9) 

n=k-M-l 


Figure  B. 2.  Number  of  sample  points  M 
for  various  passband  cutoff  widths  2  Vp 
The  upper  curve  is  for  50  dB  rejec¬ 
tion  while  the  lower  one  is  for  30  dB. 
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where  y(k)  are  filtered  samples  of  the  inphase  or  quadrature  components  x(k),and 
w(n)  are  the  weights.  This  realization  requires  2M  contiguous  input  samples  to 
produce  M  output  samples  without  transients.  Because  of  this  long  dwell  time  and 
due  to  large  storage  requirement,  the  scheme  makes  a  poor  candidate  for  interlaced 
processing. 

If  one  accepts  a  Fourier  Transform  as  a  viable  alternative  to  pulse  pair 
processing,  then  the  weights  w(k)  obtained  for  the  FIR  filter  can  be  utilized  as 
data  window  coefficient  in  which  case  the  Fourier  filter  will  have  exactly  the 
desired  characteristic. 

To  illustrate  the  point,  Hamming,  von  Hann,  and  rectangularly  windowed  time 
series  data  containing  ground  clutter  were  Fourier  analyzed  for  antenna,  rotation 
rates  of  l°s  and  10°s  .  With  a  32-point  Fourier  transform,  the  von  Hann  window 
offers  40  dB  or  more  rejection  if  three  spectral  points  (one  either  side  of  zero 
and  zeroth)  are  eliminated  (Fig.  B.3).  However,  this  type  of  filtering  implies 
a  notch  of  about  4  m'S  ,  which  may  still  be  acceptable.  The  notch  decreases  to 
about  2  m-s~  when  64  pulses  are  utilized  (Fig.  B.4). 

An  important  point  that  must  be  considered  in  using  the  weighted  Fourier 
Transform  for  clutter  rejection  is  the  desirability  to  achieve  perfect  DC 
cancellation  for  all  but  a  few  spectral  lines.  This  can  be  obtained  only  with  some 
weighting  functions  (rectangular,  von  Hann,  Hamming,  see  Figs.  B. 5  and  B. 6).  In 
general,  Fourier  Transform  of  optimum  weights  will  have  nonzero  lobes  at  the 
location  of  Fourier  coefficients  if  these  are  equally  spaced  such  as  in  the  FFT. 
This  implies  that  a  pure  DC  line  would  spill  into  other  coefficients  when  FFT  is 
used  to  calculate  the  transform.  A  Discrete  Fourier  Transform  (DFT)  can  be  chosen 

such  that  the  coefficients  are  in  the  zeroes  of  the  lobes, but  the  computations  are 
lengthy. 
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Figure  B.  3.  Power  spectra  of  ground  clutter  at  antenna  rotation  rates  of  l°s~ 
(left  graphs)  and  10°s~l  'right  graphs).  Data  windows  applied  to  time  series 
are  indicated.  Number  of  samples  is  32.  Solid  move  connects  the  data  points 
of  the  main  peak  (lobe)  while  in  the  noise  it  is  a  five  point  running  average. 
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NORMALIZED  FREQUENCY  OR  VELOCITY  (fTt,  v/2v0) 
01  23456789  10 


Figure  B.  5.  Spectral  window  corresponding  to  the  von  Hann  data  window 


NORMALIZED  FREQUENCY  OR  VELOCITY  (fT,,v/2v0) 
01  23456789  10 


Figure  B.6.  Spectral  window  corresponding  to  Homing  data  window. 
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APPENDIX  C 

Clutter  Rejection  in  the  Autocovariance  Domain 


The  ground  clutter  echoes  have  generally  narrow  Gaussian  shaped  power  spectrum 
and  are  nearly  always  centered  on  zero  Doppler  frequency.  Therefore,  it  may  be 
feasible  to  remove  the  clutter  autocorrelation  function  from  the  real  part  of  the 
composite  autocorrelation  prior  to  calculation  of  power,  velocity,  and  spectrum 
width  by  pulse  pair  algorithm.  We  refer  to  this  technique  as  the  autocovariance 
technique  for  obtaining  spectral  moment  estimates  in  the  presence  of  ground 
clutter. 

1.  Power  Spectrum  Estimation 

A 

The  total  echo  power  P  in  the  presence  of  clutter  can  be  defined  as 


p  =  S  +  C*N*jjf  $  (I^+Q,2)  'C-1) 

where  I..  and  Q.  are  the  quadrature  video  components. 

It  can  be  shown  that  for  an  assumed  Gaussian  spectrum,  the  clutter  power 
estimate  may  be  obtained  from  the  equation 


i  i. 

«  Ui  ’ 


i2  +  < 


lC.2) 


where  wc  is  the  clutter  spectrum  width  estimate  in  (Hz)  as  opposed  to  (m-s“^), 

M  is  the  number  of  sample  points;  Ts  is  the  time  lag,  and  terms  in  parenthesis 
are  the  DC  values  squared.  To  estimate  wc,  we  may  use  the  known  antenna  rotation 
rate  and  equation  (1.6).  Equation  (C.l)  was  programed  using  simulated  data  and 
estimates  of  logarithmic  signal  power  (log  §)  for  different  clutter  to  signal  ratios 
(C/S)  are  plotted  on  Fig.  C.l.  In  these,  simulations  in  lieu  of  an  estimated  wc, 
a  true  input  clutter  spectrum  width,  was  used.  At  C/S  =  0  dB  and  no  filtering, 
a  bias  error  of  3  dB  in  the  power  estimate  is  obtained  because  the  power  is  a 
sum  of  clutter  and  signal  power.  However,  with  the  partial  compensation  via  the 
autocovariance  technique  (C.l)  this  bias  is  reduced  to  2  dB.  This  is  not  a 
significant  reduction  and  thus  the  technique  cannot  be  recommended  for  power 
estimation  in  presence  of  ground  clutter. 

2.  Mean  Velocity  Estimation 

The  total  autocovariance  estimate  R.(T  )  for  a  signal  and  clutter  that  are  not 
correlated  is  the  sum  of  the  individual  autocovariances: 


«t<v  ■  Rs<v  +  MV 


(C.3) 


A  A 

where  R  (T  )  and  R  (T  )  are  the  estimated  signal  and  clut  »r  autocovariances. 

J  J  V  b 

The  pulse  pair  velocity  estimate  in  the  presence  of  tne  clutter  is  defined  as 


arg 


Im[Rs(Ts)1 

Re[Rs(Ts)] 


IC.4) 


where 


ImtRs(Tsn  =  Im  [Rt(Ts)]  -  Im  [RC<T,.)]  (C.5) 

Since  the  mean  of  R  (T  )  is  real,  the  following  relationships  hold  between  the  mean- 
or  true  values:  c  s 

Im[Rs(Ts)]  =  Im[Rt(Ts)]  IC.6) 


and 


Re[Rs(TsJ]  =  Re[Rt  ITSJ]  -  Re[Rc(Ts)]  (C.7) 

These  equalities  are  assumed  to  hold  also  for  the  respective  estimates. 

The  real  and  imaginary  part  of  the  total  estimated  autocovariance,  Rt(T$),  can 
be  represented  in  terms  of  the  inphase  and  quadrative  components  as 


RG[VTs>]  =  M  l=]^i+l  !i  +  Qi+1  Qi 5 

taC*t<Ts»-ff  !=1  Wr1!*!  V 


CC.8) 


The  autocovariance  estimate  for  Gaussian  shaped  ground  clutter,  RC(TS ) ,  is 


where,  $  is  the  clutter  spectrum  width  estimate,  (m-s-1).  Thus,  equation  (C.4) 
in  termscof  the  estimated  velocity  becomes: 


Ini[rt.  (T  / 

•*>  — 

S 


2  2  /0  2 1 
-  -  /2v,  I 

Re[Rt(Ts)]  -  C  e  c  a 


(C.10) 


The  above  equation  was  tested  on  simulated  data  and  estimates  of  velocity  for 
different  ground  clutter  to  signal  ratios  C/S  are  plotted  on  Fig.  C.2. 

When  C/S  =  0  dB,  j  =  0.2  m.s“\  and  1.0  m.s”^<  a  <  8.0  m.s"\  the  bias  of 
the  mean  velocity  estimate  decreases  by  33  percent  which  translates  to  an  effective 
S/C  improvement  of  about  3-5  dB.  For  C/S  =  10,  20,  and  30  dB,  a  heavily  biased 
estimate  of  the  nun  velocity  is  obtained  which  indicates  that  the  autocovariance 
technique  does  not  really  perform  well  when  C/S  >  0  dB.  Estimated  mean  velocities 
for  a  =  .1  m-s"  and  a  =  .3_m-s”  were  also  obtained.  The  results  are  similar 
to  thi  ones  when  a  =  0;2  m.s"  and  thus  are  not  presented  here.  The  algorithm 
did  poorly  for  thecclutter  width  values  a  -  0.3  m.s"  .  Note  that  the  pulse  Dair 
processor  generates  unbiased  mean  velocity  estimates  (if  no  filtering  is  done) 
only  for  C/S  <  -10  dB  in  agreement  with  theoretical  predictions  of  Appendix  A. 

3.  Spectrum  Width  Estimation 

For  Gaussian  spectra,  the  spectrum  width  estimate  in  the  presence  of  ground 
clutter  is  defined  as 


ttT 


iyv! 


sgr.|sn 


Rs<V 


(c.ll) 


where  S  is  the  signal  power  as  defined  in  equation  (C.ll  and  | ^  (T  )[  is  defined 
with  (C.6)  and  C.7).  Using  equations  (C.8)  and  (C.9)  | R  (T  ) |  scan  be  represented 
in  discrete  form  as  / 


MV  -  [m  f^iV  +(VV 


-  C  e 


-V^aT 


'  M  I2 

+  ^  i=i(Qi+lIi  "  w 

through  the  estimator  for  diff 


(C. 12) 


Simulated  data  were  run  through  the  estimator  for  different  ground  clutter  to 
signal  ratios  (C/S).  Figure  (C.3)  shows  the  plot  of  estimated  mean  width  versus 
true  mean  width  without  filtering  and  for  the  autocovariance  filtering. 


Comparison  of  the  two  cases  for  a  C/S  =  0  dB  reveals  that  the  bias  of  the 
spectrum  width  estimate  decreases  by  about  50  percent  when  the  autocovariance 
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technique  is  used.  But  even  then  the  bias  is  significant  and  what's  more  it  does 
depend  on  both  the  mean  velocity  and  clutter  to  signal  ratio.  This  means  the 
technique  is  unacceptable  for  spectrum  width  estimation.  Moreover,  we  see  that 
accurate  estimation  of  spectrum  widths  with  the  pulse  pair  technique  requires 
larger  S/C  ratio  (15  dB)  than  the  mean  velocity  estimation  (10  dB)  in_^greement 
with  results  of  Appendix  A.  Clutter  width  values  larger  than  0.3  m-s  were  also 
tried  but  the  algorithm  failed  likewise. 
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Figure  C.l.  Logarithm  of  estimated  mean  power  versus  true  mean  velocity  for  no 
filtering  (left  figure)  and  autocovariance  technique  (right  figure).  Clutter 
spectrum  width  CLWD  of  .2  m-s'1  is  indicated,  and  so  is  signal  width  SIGWD  = 
1.0  m-s'1.  Clutter  to  signal  ratios  range  from  0  to  -<*>  dB  in  the  no  filter 
graph  and  0  to  30  dB  ir.  the  autocovariance  graph. 
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Figure  C.2  Plot  of  estimated  mean  velocity  versus  true  mean  velocity  for  no 
filtering  (upper  graphs)  and  autocovariance  technique  (lower  graphs )• 
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Figure  C. 3  Plot  of  estimated  mean  width  versus  true  mean  width  with  no  filtering 
(top  graphs)  and  with  autocc">ariance  filtering  (bottom  graphs).  Mean  velocities 
of  the  signal  spectrum  (4,  12,  and  24  m-s~1)  as  well  as  the  clutter  to  signal 
ratios  are  indicated.  Note  that  C/S  ratios  are  less  than  0  dB  for  no  filtering 
and  are  larger  than  0  dB  for  autocovariance  filtering. 


In  summary,  simply  removing  the  ground  clutter  from  the  composite  autocovariance 
improves  very  little  the  estimates  of  reflectivity,  mean  velocity  and  Doppler  spec¬ 
trum  width.  The  improvement  is  equivalent  to  about  3-5  dB  increase  in  signal  to 
clutter  power  ratios.  Therefore,  this  technique  is  not  suitable  for  Doppler  weather 
radars . 
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APPENDIX  D 

Antenna  Speed,  Apparent  Pattern,  and  Dwell  Time 


We  have  examined  the  various  rotational  rates  that  could  apply  to  a  radar  such 
as  the  NSSL  Doppler.  The  philosophy  of  design  and  curves  depicting  the  errors  are 
listed. 

Let's  start  with  the  beamwidth,  0,  =  .8°  which  is  the  beamwidth  of  NSSL  Doppler 
radars.  Now  if  we  agree  that  the  resolution  should  not  be  less  than  1.2°  (i.e., 

50  percent  worse  than  beamwidth),  we  can  find  the  relationship  between  the  antenna 
rotation  a,  dwell  time  MT  ,  and  apparent  beamwidth  <|>  .  From  Fig.  D.l  (Zrnic1  and 
Doviak,  1976)  we  see  thats  a 


$g/®i  =  1 -5 


(D.la) 


aMye,  =  l 


(D.l b; 


Thus,  Eq.  D.lb  is  basic  as  it  ties  the  rotation  rate  and  dwell  time  to  the  beamwidth. 
As  illustration,  consider  the  following  three  rotation  rates:  18°/s  (3  rpm), 

36°/s,  and  54°/s. 

1.  Velocity  Estimation 

On  Fig.  D.2a  and  b  art  plotted  the  standard  deviation  of  velocity  estimates 
for  a  radar  transmitting  contiguous  pulses  with  parameters  as  at  NSSL.  Auto¬ 
covariance  processing  is  assumed. 

We  note  from  Fig.  D.2b  that  at  9  rpm  a  's  up  to  5.25  rn-s"1  produce  errors  less 
than  2  m-s  .  This  indicates  that  there  is  no  need  for  exotic  improvements  in  the 
signal  design  and  processing  algorithms  for  mean  velocity  estimation. 

2.  Reflectivity  Estimation 

This  subje.ct  has  been  treated  intensely  in  the  literature  and  suffices  to 
state  here  that  the  number  of  independent  samples  for  large  M  is  (Zrnic1  1979): 


Mj  =  M2tt  °vTs/(A.3-1.34) 


(D.2) 


Let  T  =  768  ys,  \  =  ,1m  and  examine  the  case  of  M=57  (i.e.,  rotation  rate  equals 
18°/s7.  Substitution  of  those  values  into  (p,2)  gives: 


Mj  =  Mo  -6.7-lU"^  =  a  *3.8 


(D.  3) 


Then  the  SD  of  reflectivity  Z  is: 


$D(Z)  =  5.57//[o" 


(D.4) 


1 


Take  a  =  1  m-s"  as  the  smallest  width  and  we  get  SD(Z)  ~  2.86  dB.  This  value 
can  bevfurther  reduced  by  averaging  in  range.  For  instance,  if  6  samples  are 
averaged,  we  obtain 

SDg(Z)  =  2.86//S-  =  1.16  dB  (0.5) 

In  order  to  have  the  SD  below  1  dB,  one  must  average  about  8  samples  in  range. 

This  illustrates  that  the  reflectivity  estimation  imposes  a  more  stringent  require¬ 
ment  to  the  number  of  samples  and  rotation  rate  than  the  velocity  estimation! 

3.  Spectrum  Width  Estimation 

Conclusions  concerning  spectrum  width  are  compatible  with  the  ones  concerning 
the  mean  velocity.  Figures  D.2c  and  d  illustrate  the  standard  deviations  for  two 
different  SNR's.  Therefore,  a  3  rpm  or  even  higher  rotation  rate  can  be  achieved 
with  little  degradation  in  the  estimates  and  with  simple  contiguous  processing 
of  samples. 


Figure  D. 1  Apparent  beamwidth  nor¬ 
malised  to  the  3  dB  bearmidth  versus 
the  normalized  rotation  rate  (aMT  /QJ. 
A  SO  percent  increase  of  the  apparent 
beamwidth  with  respect  to  the  3  dB 
beamwidth  occurs  at  vMT  /8t=2. 
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APPENDIX  E 
Program  Flow  Charts 


This  Appendix  contains  flow  charts  of  various  programs  that  were  used  throughout 
this  study.  The  most  important  subroutines  are  described  in  the  few  paragraphs 
that  follow: 


ESTPP 

This  subroutine  calculates  the  signal  power, 
velocity,  and  spectrum  width  using  pulse-pair 
algorithm. 

ESTCOV 

Signal  power,  velocity,  and  spectrum  width  are 
calculated  using  autocovariance  techniques 
(outlined  in  Appendix  C). 

IQREAD 

This  subroutine  combines  the  51  blocks  (each 
containing  5120  bytes,  4096  samples  for  the 

16  gates)  into  one  262,190  bytes  block. 

IQDECODE 

I  and  Q  time  series  data  for  each  gate  are  obtained 
from  the  262,190  bytes,  i.e.,  4096  sample  points. 

ELLIPSE 

Third  order  elliptic  filter. 

AGC 

Every  64  sample  points  are  successively  multiplied 
and  divided  by  two. 

TIMDOP 

It  generates  simulated  I  and  Q  time  series  data 
with  a  specified  mean  velocity,  CNR,  SNR,  oQ  and  ay 

DPVLWD 

From  the  input  I  and  Q  time  series  data  it  computes 
the  power,  velocity,  and  spectrum  width  using 
pulse-pair  algorithm  and  FFT. 

DPDSPY 

SHIFTIQ 

WINDOW 

SMOOTH 

IQPLOT 


Prints  the  weather  signal  power  spectrum  on  the 
printer  from  the  I  and  Q  time  series  data. 

It  eliminates  specified  number  of  time  series  data 
points  from  the  beginning  of  the  series  after  filter¬ 
ing  in  order  to  reduce  the  transient  effects. 

The  I  and  Q  time  series  data  are  weighted  with 
the  von  Hann  (raised  cosine)  or  Hamming  window. 

The  first  forty  points  of  the  simulated  I  and  Q 
time  series  data  are  smoothed  with  the  raised 
cosine  before  the  entire  8192  sample  points  are 
run  through  the  filter. 

Plots  I  and  Q  time  series  data  on  the  printer. 
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REAL  AVS(5) ,ASIGWID(4) ,AGCNR(4) ,AGCN(3) ,ZIQ(2,8192) 
NAVG  =  70  MSAM  =  128  SNR  =  20 


DO  IGCW 

_ i _ - — . — — - - 

=  1,3 

DO 

IVS  =  1,5 

DO  ISIGWID  =  1 ,4 

DO  IGCNR  =  1 ,4 

N1  =  8192 

CALL  TIMD0P(ZIQ,N1 ,VS,SIGWID,SNR,NORMAL,GCNR,GCWID) 
I  =  1 

CALL  SMOOTH(ZIQ( 1 ,1), 8192,40) 

CALL  ELLIPSE(IFILTER,ZIQ(1 ,1) ,8192) 

CALL  SHIFTIQ(ZIQ(1  «.I),ZIQ(1 ,1)  ,8064,128) 


DO  f!V  =  l.NAVG 

CALL  ESTPP(ZIQ(1 ,1) ,64,PARM(1) ,PARM(2) ,PARM(3) ,GCWID,SNR) 
1=1+  MSAM 


CALL  TIMDOP(  . 
I  =  1 

CALL  SMOOTH (  . 
CALL  ELLIPSE ( 
CALL  SHIFTIQ( 


END  DO 

PRINT  NAVG ,GCMID ,SNR , VS ,SIGWID ,GCNR ,P , V ,W , 


END  DQ 

END 

DO 

END 

DO 

END  DO 

START 


REAL  AVS(5) ,ASIGWID(4) ,AGCNR(4) ,AGCW(3) ,ZIQ(2 ,8192) 
NAVG  =  70  _ MSAM  =  128 _ SNR  =  20 

DO  IGCW  =  1 ,3 
DO  IVS  =  1 ,5 

DO  ISIGWID  =  1 ,4 
DO  IGCNR  -1,4 

N1  =  8192 

CALL  TIMD0P(ZIQ,N1  ,VS,SIGWID,SNR,NORMAL,GCNR,GCWID) 
I  =  1 

CALL  SHI FTIQ(ZIQ( 1 ,1) ,ZIQ(1 ,1) ,8064,128) 


DO  NV  =  1 ,NAVG 

CALL  ESTCOV ( ZIQ( 1 ,1) ,64,PARM(1 ) ,PARM(2) ,PARM(3) ,GCWID,SNR) 
1=1+  MSAM 


CALL  timDOP(  . 

I  =  1 

CALL  SHIFTIQ(  .  ) 


ND  DO 

PRINT  NAVG ,GCWID ,SNR ,VS ,SIGWID ,GCNR ,P ,V ,W 


END  DO 

END 

DO 

END 

DO 

END  DO 

START 


REAL  AVS(5) ,  ASIGWID(5),  AGCNR(2) ,  AGCW(l),  Ziq(2,8T92) 
MSAM  =  64  NAVG  =  100  SNR  =  20 


DO  IGCW 

=  1,1 

DO 

IVS  =  1,5 

DO  ISIGWID  =  1 ,5 

DO  IGCNR  =  1,2 

N1  =  8192 

CALL  TIMD0P  (ZIQ.N1 ,VS,$IGWID, SNR, NORMAL, GCNR,GCWID) 
I  =  1 


DO  NV  =  UNAVG 

CALL  ELLIPSE(IFILT£R,ZIQ(1 ,1) ,MSAM) 

CALL  ESTPP(ZIQ(1 ,1) ,64,PARM(1 } ,PARM(2),PARM(3) ,GCWID,SNR) 


I  =  J.  +  MSAM 


I:8192j>- 


_> _ j  CALL  TIMDOP  ( 

I  =  1 


END  DO 

PRINT  NAVG, GCWID, SNR, SIGWID,GCNR,P,V,W 


END  DO 

FND 

DO 

END 

DO 

END  DO 

.•1*1 


REAL  AVS ( 5 ) ,  ASI6WID(5),  AGCNR(2),  AGCW(1),  ZIQ(2,8192) 
AMPD  =  2  MSAM  =  64  NAVG  =  100  SNR  =  20 


DO  IGCW 

=  1,1 

DO 

IVS  =  1,5 

DO  ISIGWID  =1,5 

DO  IGCNR  =1,2 

N1  =  8192 

CALL  TIMD0P  (ZIQ.Nl , VS, SIGWID, SNR, NORMAL, GCNR.GCWID) 


I  =  1 

CALL  AGC  (ZIQ(l.I),  N1 ,  AMPD) 

CALL  ELLIPSE (I  FI LTER,ZIQ(1 ,1) ,N1 ) 


DO  NV  =  1 ,  NAVG 

CALL  ESTPP(ZI0(1 .1 ) .64,PARM(1 ) ,PARM(2) ,PARM(3) .GCWID.SNR) 


1=1+  MSAM 

1:81 


CALL  TIMDOP  ( 
I  = 


END  DO  _ 

PRINT  NAVG,GCWID:SNR,SIGWID,GCNR,P.V,W 


END  DO 
END  DO 
END  DO 
END  DO 


CHAR  RECORD  (262200) 

REAL  ZIQ(2,4Q96) ,  P0WER(128) 


READ  INITIAL  AND  ENDING  TIME 
CALL  IOREAD  (7, RECORD,  ITYPE) 

CALL  IRDECODE  ( RECORD, NSAM,NDATE,NTIME,AZM»ELEV, 
NSTEP , RNGMS , VMAX , ZI 0) 


IGATE  1  =  2 

IGATE  2  =  3 

DO  I  =  I  GATE! ,IGATE2 

CALL  IODECODE  (  ....  I) 


MSAM  =  128  SNR  =  100 

CALL  DPVLWD(ZIQ,  MSAM.  SNR.VELMAX,POWER,VPPS.WPPS.VFFT,WFFT) 
CALL  DPDSPY{ POWER, MSAM,  ' HANN * ) 

CALL  IQPLOT (ZIO,  MSAM,  1) 

CALL  ELLIPSE(ZIQ,  NSAM) 


K  =  NSAM  -  81  1 

CALL 

SHIFTJQ(ZIQ.  ZIQ,  K.ISHIFT) 

CALL 

DPVLWD(ZIQ,MSAM, . . .) 

CALL 

DPDSPY  (  .  ) 

CALL 

IQPLOT  (  .  ) 

END 

DC 
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REAL  ZIQ(2,256),  XIQ(2,123) 
200  READ  IBTM. IETM, DELAY! ,  DELAY 2 


1GATE  1  =  2 
IGATE  2=3 

DO  300  I  =  IGATE 1 ,  IGATE2 
MSAM  =  128  SNR  =  100 


CALL  MO VEARRY ( Z I Q , X I 0 , MSAM ) 

A  =  B  =  0.5 

CALL  WINDOW(ZIO,MSAM,A,B) 

CALL  DPVLWD(ZIQ, MSAM, SNR, VELMAX, POWER, VPPS,WPPS,VFFT,WFFT) 
CALL  DPDSPY (POWER,  MSAM.  1 HANN ' ) 

CALL  IQPLOTfZIQ,  MSAM,  1) 


CALL 

MOVEARRY ( X I Q , Z I 0 , MSAM ) 

CALL 

SMOOTH ( Z I Q , MSAM , 40 ) 

CALL 

ELLIPSE(ZIQ,MSAM) 

CALL 

SHIFTIQ(ZI0,ZIQ,64,64) 

MSAM 

-  64 

A  =  B 

=  0.5 

CALL 

WINDOW(ZIO,MSAM,A,B) 

CALL 

DPVLWD(ZIQ,MSAM . ) 

CALL 

DPDSPY (POWER, MSAM,  'HANN') 

CALL 

IQPLOT(ZIQ.MSAM,  1) 
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Two  versions  of  the  elliptic  filter  subroutine  are  listed.  The  first  one  uses 
initialization  to  anticipate  steady  statt.  The  second  one  initializes  memory 
elements  to  zero.  The  first  subscript  (i.e.,  1)  refers  to  I  while  2  refers  to  the 
0  coripcne'u.  Z's  are  input  time  samples;  y's  are  filtered  time  samples.  XK's 
are  filter  coefficients  corresponding  to  Fig.  3.1  and  Table  3.1. 

^  Subroutine  ELLIPSE  1  (IFILTER^.N) 

REAL  Z(2.N).y(2,819?) *XK(4) 


__.ij-LAiL'iinvr  nriTri,ivnf^t“-*-J^‘A‘-  • 


ST 

I 


I 
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